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Abstract. Recent experimental and theoretical studies show that the switching 
behavior of magnetic nanoparticles can be well controlled by external time-dependent 
magnetic fields. In this work, we inspect theoretically the influence of the temperature 
and the magnetic anisotropy on the spin-dynamics and the switching properties of 
single domain magnetic nanoparticles (Stoner-particles). Our theoretical tools are 
the Landau-Lifshitz-Gilbert equation extended as to deal with finite temperatures 
within a Langevine framework. Physical quantities of interest are the minimum 
field amplitudes required for switching and the corresponding reversal times of the 
nanoparticle's magnetic moment. In particular, we contrast the cases of static and 
time-dependent external fields and analyze the influence of damping for a uniaxial and 
a cubic anisotropy. 
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1. Introduction 

In recent years, there has been a surge of research activities focused on the spin dynamics 
and the switching behavior of magnetic nanoparticles pQ. These studies are driven 
by potential applications in mass-storage media and fast magneto-electronic devices. 
In principle, various techniques are currently available for controlling or reversing the 
magnetization of a nanoparticle. To name but a few, the magnetization can be reversed 
by a short laser pulse [2], a spin-polarized electric current (3j H] or an alternating 
magnetic field [HI El d [H [91 [THl HH H21 [13]. Recently ©, it has been shown for a 
uniaxial anisotropy that the utilization of a weak time- dependent magnetic field achieves 
a magnetization reversal faster than in the case of a static magnetic field. For this case 
[6], however, the influence of the temperature and the different types of anisotropy 
on the various dependencies of the reversal process have not been addressed. These 
issues, which are the topic of this present work, are of great importance since, e.g. 
thermal activation affects decisively the stability of the magnetization, in particular 
when approaching the superparamagnetic limit, which restricts the density of data 
storage [14J. Here we study the possibility of fast switching at finite temperature with 
weak external fields. We consider magnetic nanoparticles with an appropriate size as to 
display a long-range magnetic order and to be in a single domain remanent state (Stoner- 
particles). Uniaxial and cubic anisotropies are considered and shown to decisively 
influence the switching dynamics. Numerical results are presented and analyzed for 
iron-platinum nanoparticles. In principle, the inclusion of finite temperatures in spin- 
dynamics studies is well-established (cf. (191 EDI [23], [151 EE [I] an d references therein) 
and will be followed here by treating finite temperatures on the level of Langevine 
dynamics. For the analysis of switching behaviour the Stoner and Wohlfarth model 
(SW) [T7] is often employed. SW investigated the energetically metastable and stable 
position of the magnetization of a single domain particle with uniaxial anisotropy in 
the presence of an external magnetic field. They showed that the minimum static 
magnetic field (generally referred to as the Stoner- Wohlfarth (SW) field or limit) needed 
to coherently reverse the magnetization is dependent on the direction of the applied 
field with respect to the easy axis. This dependence is described by the so-called 
Stoner- Wohlfarth astroid. The SW findings rely, however, on a static model at zero 
temperature. Application of a time- dependent magnetic field reduces the required 
minimum switching field amplitude below the SW limit [6j. It was, however, not yet 
clear how finite temperatures will affect these findings. To clarify this point, we utilize 
an extension of the Landau-Lifshitz-Gilbert equation [18] including finite temperatures 
on the level of Langevine dynamics [191 EDI E3]. Our analysis shows the reversal time 
to be strongly dependent on the damping, the temperature and the type of anisotropy. 
These dependencies are also exhibited to a lesser extent by the critical reversal fields. 
The paper is organized as follows: next section 2 presents details of the numerical scheme 
and the notations whereas section 3 shows numerical results and analysis for FesoPtso 
and Fe7oPt 30 nanoparticles. We then conclude with a brief summary. 
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2. Theoretical model 

In what follows we focus on systems with large spins such that their magnetic dynamics 
can be described by the classical motion of a unit vector S directed along the particle's 
magnetization n, i.e. S = fj,/ ' us and us is the particle's magnetic moment at saturation. 
The energetics of the system is given by 

n=n A + n F . (i) 

where 1~La {Hf) stands for the anisotropy (Zeeman energy) contribution. Furthermore, 
the anisotropy contribution is expressed as H,a = —Df(S) with D being the anisotropy 
constant. Explicit form of f(S) is provided below. The magnetization dynamics, i.e. the 
equation of motion for S, is governed by the Landau-Lifshitz-Gilbert (LLG) equation 



B e (t) + a(S x B e (t)) . (2) 



dt (1 + a 2 ) 

Here we introduced the effective field B e (t) = —l/(/j,s)&H/dS which contains the 
external magnetic field and the maximum anisotropy field for the uniaxial anisotropy 
Ba = IB I [is- 7 is the gyromagnetic ratio and a is the Gilbert damping parameter. The 
temperature fluctuations will be described on the level of the Langevine dynamics [19|. 
This means, a time-dependent thermal noise adds to the effective field B e (i) [19|. 
C(t) is a Gaussian distributed white noise with zero mean and vanishing time correlator 

(m<i(t)) = ^^wt-n (3) 

H'sl 

i,j are Cartesian components, T is the temperature and ks is the Boltzmann constant. 
It is convenient to express the LLG in the reduced units 

b = -f^, t = u a t, u a = jB A . (4) 

DA 

The LLG equation reads then 
dS 1 



-S x 



b(r) + a(S x b(r)) , (5) 



Or (1 + a 2 ) 
where the effective field is now given explicitly by 

with 

(Qiir^ir)) = eS id 6(r - r'); e = 2 -^f. (7) 

The reduced units are independent of the damping parameter a. In the following sections 
we use extensively the parameter 

? = -0". (8) 

q is a measure for the thermal energy in terms of the anisotropy energy. And 
d = B/(h s Ba) expresses the anisotropy constant in units of a maximum anisotropy 
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energy for the uniaxial anisotropy and is always 1/2. The stochastic LLG equation (JS]) 
in reduced units (j3J) is solved numerically using the Heun method which converges in 
quadratic mean to the solution of the LLG equation when interpreted in the sense of 
Stratonovich [20]. For each type of anisotropy we choose the time step At to be one 
thousandth part of the corresponding period of oscillations. The values of the time 
interval in not reduced units for uniaxial and cubic anisotropies are At ua = 4.61 ■ 10~ 15 s 
and At ca = 64.90 • 10~ 15 s, respectively, providing us thus with correlation times on the 
femtosecond time scale. The reason for the choice of such small time intervals is given in 
[j~9] , where it is argued that the spectrum of thermal-agitation forces may be considered 
as white up to a frequency of order ksl '/h with h being the Planck constant. This value 
corresponds to 10~ 13 s for room temperature. The total scale of time is limited by a 
thousand of such periods. Hence, we deal with around one million iteration steps for 
a switching process. Details of realization of this numerical scheme could be found in 
references [2TJ [22j [20] . We note by passing, that attempts have been made to obtain, 
under certain limitations, analytical results for finite-temperature spin dynamics using 
the Fokker-Planck equation (cf. [T5], [T6] and references therein). For the general case 
discussed here one has however to resort to fully numerical approaches. 

3. Results and interpretations 

We consider a magnetic nanoparticle in a single domain remanent state (Stoner-particle) 
with an effective anisotropy whose origin can be magnetocrystalline, magnetoelastic and 
surface anisotropy. We assume the nanoparticle to have a spherical form, neglecting 
thus the shape anisotropy contributions. In the absence of external fields, thermal 
fluctuations may still drive the system out of equilibrium. Hence, the stability of 
the system as the temperature increases becomes an important issue. The time t at 
which the magnetization of the system overcomes the energy barrier due to the thermal 
activation, also called the escape time, is given by the Arrhenius law 



where the exponent is the ratio of the anisotropy to the thermal energy. The coefficient 
to may be inferred when D ^> ksT and for high damping |19j (see [25] for a critical 
discussion) 



Here we focus on two different types of iron-platinum-nanoparticles: The compound 
FesoPtso which has a uniaxial anisotropy [26], [27], whereas the system Fe7oPt3o possesses 
a cubic anisotropy [24J. Furthermore, the temperature dependence will be studied by 
varying q (cf. eq.(JED)- 

For FesoPtso the important parameters for simulations are the diameter of the 
nanoparticles 6.3 nm, the strength of the anisotropy K u = 6 • 10 6 J/m 3 , the magnetic 
moment per particle \i v = 21518 • /i# and the Curie-temperature T c = 710K [26 :> 127) . 



D 



t = t e k B T , 



(9) 




(10) 



Temperature dependent magnetization dynamics of magnetic nanoparticles 



5 



The relation between K u and D u is D u = K U V U , where V u is the volume of FesoPtso 
nanoparticles. In the calculations for FesoPtso nanoparticles the following q values 
were chosen: q\ = 0.001, q<i = 0.005 or q% = 0.01 which correspond to the real 
temperatures 56K, 280K or 560K, respectively (these temperatures are below the 
blocking temperature). The corresponding escape times are t qi ~ 2 • 10 217 s, t Q2 « 10 75 s 
and t qs ~ 7- 10 31 s, respectively. In some cases we also show the results for an additional 
temperature g i = 0.0001 with the corresponding real temperature to be equal to 5K. 
The corresponding escape time for this is t qoi m 10 4300 s. These times should be compared 
with the measurement period which is about t m ~ 5 ns, endorsing thus the stability of 
the system during the measurements. 

For Fe7oPt3o the parameters are as follows: The diameter of the nanoparticles 2.3 nm, 
the strength of the anisotropy K c = 8 • 10 5 J/m 3 , the magnetic moment per particle 
fjL p = 2000 ■ fiB, the Curie-temperature is T c = 420K [21], and D c = K C V C (V c is the 
volume.) For Fe7oPt3o nanoparticles the values of q we choose in the simulations are 
q4 = 0.01, q$ = 0.03 or g 6 = 0.06 which means that the temperature is respectively 
0.3if , 0.9K or 1.9K. The escape times are t QA w 10 34 s, t Q5 w 2 • 10 5 s and t q6 w 2 ■ 10~ 2 s, 
respectively. Here we also choose an intermediate value go4 — 0.001 and the real 
temperature 0.03K with the corresponding escape time to be equal to t qo4 ~ 10 430 s. 
The measurement period is the same, namely about 5 ns. All values of the escape times 
were given for a = 0.1. 

Central to this study are two issues: The critical magnetic field and the corresponding 
reversal time. The critical magnetic field we define as the minimum field amplitude 
needed to completely reverse the magnetization. The reversal time is the corresponding 
time for this process. In contrast, in other studies [6] the reversal time is defined as the 
time needed for the magnetization to switch from the initial position to the position 
S x — 0, our reversal time is the time at which the magnetization reaches the very 
proximity of the antiparallel state (Fig. [T]). The difference in the definition is in so far 
important as the magnetization position S z = at finite temperatures is not stabile so 
it may switch back to the initial state due to thermal fluctuations and hence the target 
state is never reached. 

3.1. Nanoparticles having uniaxial anisotropy: Fe^Pt^ 

A FesoPtso magnetic nanoparticle has a uniaxial anisotropy whose direction defines the 
z direction. The magnetization direction S is specified by the azimuthal angle <fi and 
the polar angle 9 with respect to z. In the presence of an external field b applied at 
an arbitrarily chosen direction, the energy of the system in dimensionless units derives 
from 

H = -dcos 2 6 - S -b. (11) 

The initial state of the magnetization is chosen to be close to S z = +1 and we aim at 
the target state S z = —1. 
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Figure 1. (Color online) Magnetization reversal of a nanoparticle when a static field 
is applied at zero Kelvin (qq = 0, black) and at reduced temperature = 0.01 = 560K 
(blue). The strengths of the fields in the dimensionless units (j3| and ([5]) are b — 1.01 
and b — 0.74, respectively. The damping parameter is a — 0.1. The start position of 
the magnetization is given by the initial angle 9q — 7r/360 between the easy axis and 
the magnetization vector. 



3.1.1. Static field For an external static magnetic field applied antiparallel to the z 
direction (6 = —be z ) eq. lfTTT) becomes 

H = -dcos 2 6 + bcos9. (12) 

To determine the critical field magnitude needed for the magnetization reversal we 
proceed as follows (cf. Fig. CD): At first, the external field is increased in small steps. 
When the magnetization reversal is achieved the corresponding values of the critical 
field versus the damping parameter a are plotted as shown in the inset of Fig. [3] The 
reversal times corresponding to the critical static field amplitudes of Fig. [3] are plotted 
versus damping in Fig. HI 

In the Stoner-Wohlfarth (static) model the mechanism of magnetization reversal is not 
due to damping. It is rather caused by a change of the energy profile in the presence of 
the field. The curves displayed on the energy surface in Fig. [2] mark the magnetization 
motion in the E(6,(j)) landscape. The magnetization initiates from <f) = and 6q an d 
ends up at 9 = it. As clearly can be seen from the figure, reversal is only possible if the 
initial state is energetically higher than the target state. This "low damping" reversal 
is, however, quite slow, which will be quantified more below. For the reversal at T = 0, 
the SW-model predicts a minimum static field strength, namely b cr = B/Ba = 1 (the 
dashed line in Fig. [3]). 

This minimum field measured with respect to the anisotropy field strength does not 
depend on the damping parameter a, provided the measuring time is infinite. For T > 
the simulations were averaged over 500 cycles with the result shown in Fig. [3j The one- 
cycle data are shown in the inset. Fig. [3] evidences that with increasing temperature 
thermal fluctuations assist a weak magnetic field as to reverse the magnetization. 
Furthermore, the required critical field is increased slightly at very large and strongly 
at very small damping with the minimum critical field being at a ~ 1.0. The reason for 
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Figure 2. (Color online) The trajectories of the magnetization unit vector 
parameterized by the angles 9 and <j> at zero temperature. Other parameters are as in 
Fig. rjjfor q . 
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Figure 3. (Color online) Critical static field amplitudes vs. the damping parameters 
for different temperatures averaged over 500 times. Inset shows not averaged data for 
q 3 = 0.01 ee 560K. 



this behavior is that for low damping the second term of equation ([2]) is much smaller 
than the first one, meaning that the system exhibits a weak relaxation. In the absence 
of damping, higher fields are necessary to switch the magnetization. For high a, both 
terms in equation ([2]) become small (compared to a low-damping case) leading to a 
stiff magnetization and hence higher fields are needed to drive the magnetization. For 
moderate damping, we observe a minimum of switching fields which is due to an optimal 
interplay between precessional and damping terms. Obviously, finite temperatures do 
not influence this general trend. 

For the case of go — 0, the Landau-Lifshitz-Gilbert equation of motion can be solved 
analytically in spherical coordinates. The details of the solution can be found in Ref. 
[20J (eq. (Al)-(A8)). The final result of the solution in this reference differs, however, 
from the one given here due to to different geometries in these systems. In contrast to 
our alignment of the magnetization and the external field, the static field in Ref. [20] is 
applied parallel to the initial position of the magnetization. For the solution, we assume 
that the magnetization starts at 9 = 9q = 7r/360 and arrives at 9 = 7r. Note, that the 
expression 9 ^ is important only for zero Kelvin since the switching is not possible if 
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the magnetization starts at 9q — (the vector product in equation vanishes). The 
reversal time in the SW-limit is then given by 

1 + a 2 

Uev = g(0 , b) , (13) 

a 

where g is defined as 

From this relation we infer that switching is possible only if the applied field is larger than 
the anisotropy field and the reversal time decreases with increasing b. This conclusion is 
independent of the Stoner-Wohlfarth model and follows directly from the solution of the 
LLG equation. An illustration is shown by the dashed curve in Fig. HI which was a test 
to compare the appropriate numerical results with the analytical one. As our aim is the 
study of the reversal-time dependence on the magnetic moment and on the anisotropy 
constant, we deem the logarithmic dependence in Eq. (jl4p to be weak and write 

9(b,»s,D)*- B2 ^_ 4Dr (15) 

This relation indicates that an increase in the magnetic moment results in a decrease 
of the reversal time. The magnetic moment enters in the Zeeman energy and therefore 
the increase in magnetic moment is very similar to an increase in the magnetic field. 
An increase of the reversal time with the increasing anisotropy originates from the fact 
that the anisotropy constant determines the height of the potential barrier. Hence, the 
higher the barrier, the longer it takes for the magnetization to overcome it. 
For the other temperatures the corresponding reversal times (also averaged over 500 
cycles) are shown in Fig. HI In contrast to the case T = 0, where an appreciable 
dependence on damping is observed, the reversal times for finite temperatures show 
a weaker dependence on damping. If a — > only the precessional motion of the 
magnetization is possible and therefore t rev — > oo. At high damping the system relaxes 
on a time scale that is much shorter than the precession time, giving thus rise to an 
increase in switching times. Additionally, one can clearly observe the increase of the 
reversal times with increasing temperatures, even though these time remain on the 
nanoseconds time scale. 



3.1.2. Alternating field As was shown in Ref. P, [7J [15] theoretically and in Ref. [5] 
experimentally, a rotating alternating field with no static field being applied can also 
be used for the magnetization reversal. A circular polarized microwave field is applied 
perpendicularly to the anisotropy axis. Thus, the Hamiltonian might be written in form 
of equation (fTTjl and the applied field is 

b(t) = bo cos ute x + bo sin ute y , (16) 

where bo is the alternating field amplitude and uj is its frequency. For a switching of 
the magnetization the appropriate frequency of the applied alternating field should be 
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Damping 



Figure 4. (Color online) Reversal times corresponding to the critical static fields in 
Fig. [3] vs. damping averaged over 500 cycles. Inset shows the as-calculated numerical 
results for qs = 0.01 = 560K (one cycle). 
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Figure 5. (Color online) Magnetization reversal in a nanoparticlc using a time 
dependent field for a — 0.1 and at a zero temperature. The field strength and frequency 
in the units ((4]) are respectively &o = 0.18 and lu — w /1.93. Inset shows for this case 
the magnetization reversal for the temperature (73 = 0.01 = 560-ftT with &o = 0.17 and 
the same frequency. 



chosen. In Ref. [15] analytically and in [6] numerically a detailed analysis of the optimal 
frequency is given which is close to the precessional frequency of the system. The role 
of temperature and different types of anisotropy have not yet been addressed, to our 
knowledge. 

Fig. [5] shows our calculations for the reversal process at two different temperatures. 
In contrast to the static case, the reversal proceeds through many oscillations on a 
time scale of approximately ten picoseconds. Increasing the temperature results in an 
increase of the reversal time. 

Fig. [6] shows the trajectory of the magnetization in the E(#,0) space related to the 
case of the alternating field application. Compared with the situation depicted in Fig. 
|2l the trajectory reveals a quite delicate motion of the magnetization. It is furthermore, 
noteworthy that the alternating field amplitudes needed for the reversal (cf. Fig. [7j) are 
substantially lower than their static counterpart, meaning that the energy profile of the 
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Trajectory in Hie e spate 




Figure 6. (Color online) Trajectories followed by magnetization as specified by 9 and 
(f> for qo = 0. Other parameters are bo = 0.18, a = 0.1 and li — w a /1.93. Energy-profile 
variations due to the oscillating external field are not visible on this scale. 



system is not completely altered by the external field. 

Fig. [7J inspects the dependence of the minimum switching field amplitude on 
damping. The critical fields are obtained upon averaging over 500 cycles. The SW- 
limit lies by 1 on this scale. In contrast to the static case, the critical fields increase 
with increasing a. In the low damping regime the critical field is smaller than in the 
case of a static field. This behavior can be explained qualitatively by a resonant energy- 
absorption mechanism when the frequencies of the applied field matches the frequency of 
the system. Obviously, at very low frequencies (compared to the precessional frequency) 
the dynamics resembles the static case. 

The influence of the temperature on the minimum alternating field amplitudes is 
depicted in Fig. [7J With increasing temperatures, the minimum amplitudes become 
smaller due to an additional thermal energy pumped from the environment. The curves 
in this figure can be approached with two linear dependencies with different slopes for 
approximately a < 1 and for a > 1; for high damping it is linearly dependent on a, 
more specifically it can be shown that for high damping the critical fields behave as 



The proportionality coefficient contains the frequency of the alternating field and the 
critical angle 9. The solution (ll?p follows from the LLG equation solved for the case 
when the phase of the external field follows temporally that of the magnetization, which 
we checked numerically to be valid. 

The reversal times associated with the critical switching fields are shown in (Fig. 
E}. Qualitatively, we observe the same behavior as for the case of a static field. The 
values of the reversal times for T = are, however, significantly smaller than for the 
static case. For the same reason as in the static field case, an increased temperature 
results in an increase of the switching times. 
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Figure 7. (Color online) Critical alternating field amplitudes vs. damping for 
different temperatures averaged over 500 times. Inset shows not averaged data for 
q 3 = 0.01 = 560K. 




Figure 8. (Color online) The damping dependence of the reversal times corresponding 
to the critical field amplitudes of Fig. [7] for different temperatures. Inset shows the 
case of zero Kelvin. 



3.2. Nanoparticles with cubic anisotropy: Fe^o Pt^o 

Now we focus on another type of the anisotropy, namely a cubic anisotropy which is 
supposed to be present for Fe7 Pt 30 nanoparticles [21]. The energetics of the system is 
then described by the functional form 

Iri = —d(S%Sy + SyS^ + S^Sg) — S • b, (18) 

or in spherical coordinates 

U = -d(cos 2 sin 2 sin 4 9 + cos 2 9 sin 2 9) - S ■ b. (19) 

In contrast to the previous section, there are more local minima or in other words more 
stable states of the magnetization in the energy profile for the FeyoPt3o nanoparticles. 
It can be shown that the minimum barrier that has to be overcome is d/12 which is 
twelve times smaller than that in the case of a uniaxial anisotropy. The maximal one is 
only d/3. 

The magnetization of these nanoparticles is first relaxed to the initial state close to 
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0o = tt/4 and 9 = arccos(l/\/3), whereas in the target state it is aligned antiparallel to 
the initial one, i. e. e = 3tc/A and 9 e = tt — arccos(l/v / 3)- In order to be close to the 
starting state for the uniaxial anisotropy case we choose O — 0.2499 • n, 9$ = 0.3042 • ir. 

3.2.1. Static driving field A static field is applied antiparallel to the initial state of the 
magnetization, i.e. 

b= -b/V3(e x + e y + e z ). (20) 

In Fig. Othe trajectory of the magnetization in case of an applied static field is shown. 
Similar to the previous section the energy of the initial state lies higher than that of 
the target state. The magnetization rolls down the energy landscape to eventually 
end up by the target state. The trajectory the magnetization follows is completely 
different from the one for the uniaxial anisotropy. Fig. [10] supplements this scenario 
of the magnetization reversal by showing the time evolution of the S z vector. Because 
of the different anisotropy type, the trajectory is markedly different from the case of 
the uniaxial anisotropy and a static field. Here we show only the S z magnetization 
component even though the other components also have to be taken into account in 
order to avoid a wrong target state. 

The procedure to determine the critical field amplitudes is similar to that described in 
the previous section. In Fig. [TT] the critical fields versus the damping parameter for 
different temperatures are shown. For qo, the critical field strength is smaller than 1. 
This is consistent insofar as the maximum effective field for a cubic anisotropy is 
In principle, the critical field turns out to be constant for all a but for an infinitely large 
measuring time. Since we set this time to be about 5 nanoseconds, the critical fields 
increase for small and high damping. On the other hand, at lower temperatures smaller 
critical fields are sufficient for the (thermal activation-assisted) reversal process. 

The behaviour of the corresponding switching times presented in Fig. [12] only 
supplements the fact of too low measuring time, which is chosen as 5 ns for a better 
comparison of these results with ones for uniaxial anisotropy. Indeed, constant jumps 
in the reversal times for T = K as a function of damping can be observed. The reason 
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Time, [ns] 

Figure 10. (Color online) Magnetization reversal of a nanoparticle when a static field 
b = 0.82 is applied and for a = 0.1 at zero temperature (black). The magnetization 
reversal for a = 0.1, b = 0.22 and = 0.06 = 1.9K is shown with blue color. 
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Figure 11. (Color online) Critical static field amplitudes vs. the damping parameters 
for different temperatures averaged over 500 times. Inset shows not averaged data for 
q 6 = 0.06= 1.9K. 
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Figure 12. (Color online) Reversal times corresponding to the critical static fields of 
Fig. [TT]vs. damping averaged over 500 times. 



why the reversal times for finite temperatures are lower is as follows: The initial state 
for T = K is chosen to be very close to equilibrium. This does not happen for finite 
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Trajectory in the e «■) space 




Figure 13. (Color online) Trajectories of the magnetization vector specified by the angles 9 and <j> at 
zero temperature. The chosen parameters are bo — 0.055 and u> — tJ a /1.93, where cJ a — 2/3a; Q . 



temperatures, where the system due to thermal activation jumps out of equilibrium (cf. 
see Fig. [TO]) . 

3.2.2. Time- dependent external field Here we consider the case of an alternating field 
that rotates in the plain perpendicularly to the initial state of the magnetization. It is 
possible to switch the magnetization with a field rotating in the xy— plane but the field 
amplitudes turn out to be larger than those when the field rotates perpendicular to the 
initial state. For the energy this means that the field entering equation ( fl9l) reads 

b(t) = (bo cos LOit cos (fio + b sin u)\t sin </> cos 9o)e x 

+ (—60 cosc^it sin 0o + bo sincuit cos0o cos^e^ + (—bo sin 9q sinu;it)e z , (21) 

where b is the alternating field amplitude and oj\ is the frequency associated with the 
field. This expression is derived upon a rotation of the field plane by the angles <po = tt/4 
and 0q = arccos(l/v / 3)- 

The magnetization trajectories depicted in Fig. [13] reveal two interesting features: 
Firstly, particularly for small damping, the energy profile changes very slightly (due 
to the smallness of bo) while energy is pumped into the system during many cycles. 
Secondly, the system switches mostly in the vicinity of local minima to acquire eventually 
the target state. Fig. [TU hints on the complex character of the magnetization dynamics 
in this case. As in the static field case with a cubic anisotropy the critical field 
amplitudes shown in Fig. [151 are smaller than those for a uniaxial anisotropy. Obviously, 
the reason is that the potential barrier associated with this anisotropy is smaller in this 
case, giving rise to smaller amplitudes. As before an increase in temperature leads to a 
decrease in the critical fields. 

The reversal times shown in Fig. [T6l exhibit the same feature as in the cases for uniaxial 
anisotropy: With increasing temperatures the corresponding reversal times increase. A 
physically convincing explanation of the (numerically stable) oscillations for the reversal 
times is still outstanding. 
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Figure 14. (Color online) Magnetization reversal in a nanoparticle using a time dependent field for 
a = 0.1 and go (black) and for q$ = 0.06 = 1.9K (blue). Other parameters are as in Fig. [T3l 
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Figure 15. (Color online) Critical alternating field amplitudes vs. damping for different temperatures 
averaged over 500 cycles. Inset shows the single cycle data at q 6 = 0.06 = 1.9K. 




Damping 



Figure 16. (Color online) The damping dependence of the reversal times corresponding to the critical 
fields of the Fig. [15] for different temperatures averaged over 500 runs. Inset shows the T = case. 

4. Summary 

In this work we studied the critical field amplitudes required for the magnetization 
switching of Stoner nanoparticles and derived the corresponding reversal times for 
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static and alternating fields for two different types of anisotropies. The general trends 
for all examples discussed here can be summarized as follows: Firstly, increasing the 
temperature results in a decrease of all critical fields regardless of the anisotropy type. 
Anisotropy effects decline with increasing temperatures making it easier to switch the 
magnetization. Secondly, elevating the temperature increases the corresponding reversal 
times. Thirdly, the same trends are observed for different temperatures: The critical field 
amplitudes for a static field depend only slightly on a, whereas the critical alternating 
field amplitudes exhibit a pronounced dependence on damping. In the case of a uniaxial 
anisotropy we find the critical alternating field amplitudes to be smaller than those for a 
static field, especially in the low damping regime and for finite temperatures. Compared 
with a static field, alternating fields lead to smaller switching times (T = K). However, 
this is not the case for the cubic anisotropy. The markedly different trajectories for the 
two kinds of anisotropies endorse the qualitatively different magnetization dynamics. 
In particular, one may see that for a cubic anisotropy and for an alternating field 
the magnetization reversal takes place through the local minima leading to smaller 
amplitudes of the applied field. Generally, a cubic anisotropy is smaller than the uniaxial 
one giving rise to smaller slope of critical fields, i.e. smaller alternating field amplitudes. 
It is useful to contrast our results with those of Ref. [15]. Our reversal times for 
AC-fields increase with increasing temperatures. This is not in contradiction with the 
findings of [15] insofar as we calculate the switching fields at first, and then deduce the 
corresponding reversal times. If the switching fields are kept constant while increasing 
the temperature [15] the corresponding reversal times decrease. We note here that 
experimentally known values of the damping parameter are, to our knowledge, not 
larger than 0.2. The reason why we go beyond this value is twofold. Firstly, the values 
of damping are only well known for thin ferromagnetic films and it is not clear how to 
extend them to magnetic nanoparticles. For instance, in FMR experiments damping 
values are obtained from the widths of the corresponding curves of absorption. The 
curves for nanoparticles can be broader due to randomly oriented easy anisotropy axes 
and, hence, the values of damping could be larger than they actually are. Secondly, due 
to a very strong dependence of the critical AC-fields (Fig. [7J e.g.) they can even be larger 
than static field amplitudes. This makes the time-dependent field disadvantageous for 
switching in an extreme high damping regime. 

Finally, as can be seen from all simulations, the corresponding reversal times are much 
more sensitive a quantity than their critical fields. This follows from the expression ( TL3l) . 
where a slight change in the magnetic field b leads to a sizable difference in the reversal 
time. This circumstance is the basis for our choice to average all the reversal times and 
fields over many times. This is also desirable in view of an experimental realization, for 
example, in FMR experiments or using a SQUID technique quantities like critical fields 
and their reversal times are averaged over thousands of times. The results presented 
in this paper are of relevance to the heat-assisted magnetic recording, e.g. using a 
laser source. Our calculations do not specify the source of thermal excitations but 
they capture the spin dynamics and switching behaviour of the system upon thermal 
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excitations. 
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